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Abstract: Rooftop distributed photovoltaic projects have been quickly proposed in China because 
of policy promotion. Before, the rooftops of the shopping mall had not been occupied, and it was 
urged to have a decision-making framework to select suitable shopping mall photovoltaic plans. 
However, a traditional multi-criteria decision-making (MCDM) method failed to solve this issue at 
the same time, due to the following three defects: the interactions problems between the criteria, 
the loss of evaluation information in the conversion process, and the compensation problems between 
diverse criteria. In this paper, an integrated MCDM framework was proposed to address these 
problems. First of all, the compositive evaluation index was constructed, and the application of 
decision-making trial and evaluation laboratory (DEMATEL) method helped analyze the internal 
influence and connection behind each criterion. Then, the interval-valued neutrosophic set was 
utilized to express the imperfect knowledge of experts group and avoid the information loss. Next, 
an extended elimination et choice translation reality (ELECTRE) II] method was applied, and it 
succeed in avoiding the compensation problem and obtaining the scientific result. The integrated 
method used maintained symmetry in the solar photovoltaic (PV) investment. Last but not least, 
a comparative analysis using Technique for Order Preference by Similarity to an Ideal Solution 
(TOPSIS) method and VIKOR method was carried out, and alternative plan X1 ranks first at the same. 
The outcome certified the correctness and rationality of the results obtained in this study. 


Keywords: shopping mall; photovoltaic plan; decision-making trial and evaluation laboratory 
(DEMATEL); interval-valued neutrosophic set; extended ELECTRE III; symmetry 


1. Introduction 


The frequent occurrence of fog or haze and other negative types of climate change in recent 
decades is the grave reality that the whole world is experiencing. The pivotal reason behind these 
environmental problems is atmospheric pollutants and greenhouse gas emissions, mainly produced by 
fossil fuel consumption. Fossil fuel supplies approximately 80% of the world’s energy, and it is drying 
up with the rapid increase of the world energy demand [1]. 

To face this situation, many countries have endorsed policies to submit fossil fuel utilization 
with renewable energy generation. Among diverse types of alternative energy, solar photovoltaic 
(PV hereinafter) energy is recognized as promising, since sunlight is unlimited and widespread and 
the converting efficiencies of photovoltaic are getting higher and higher while the manufacturing costs 
are becoming lower and lower [2]. 
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The past five years has witnessed the astonishing increase in installed cumulative globe solar PV 
capacity, which nearly quintupled from 70 GW in 2011, to 275 GW in 2016. China is following the 
worldwide trend with solar PV rapid development and gained the number one cumulative PV capacity 
of 65.57 GW in 2016. It is worth noting that the large-scale ground PV power station, newly installed, 
had 28.45 GW capacity this year, and grew by 75% compared to the last year, accounting for 89% of 
all new PV power plants in 2016, while the distributed PV, newly installed, with 3.66 GW capacity 
this year, increased by 45% in comparison to last year, and accounted for 11% of capacity of new 
installations in 2016. The scale of distributed PV development is significantly lower than that of the 
large-scale ground PV. 

Under these circumstances, the China authorities have launched the feed in tariff adjustment. 
The feed in tariff of ground solar PV generation will decrease to some extent, but in contrast to that 
of distributed PV, will not decrease at all. As shown in Figure 1, according to the 13th five-year 
(2016-2020) solar energy planning objectives of China, the goal is to build a total installed capacity 
of 150 GW solar PV, in which more than 40% of the new installed capacity will be from distributed 
PV, and to build 100 distributed PV demonstration areas. In fact, China is a country with high 
potential of solar radiation, and of generous policy subsidies to promote achieving the ambitious 
plan. Obviously, in China, distributed solar PV generation is government encouraged, well-resourced, 
environmentally friendly, and closely following the world trend. It is worth considering investment in 
such a promising project. 
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Figure 1. Chinese photovoltaic power structure, 2012-2016. 


Numerous vacant roofs of building or structures in the cities provide the best-fitting location 
for distributed solar PV, and farsighted investors have been preempting outstanding roof resources. 
There are various types of buildings, such as government building, hospitals, schools, coliseums, 
or residential and industrial buildings. Among all these types, the shopping centers possesses plenty 
of advantages, superior to the other types of buildings, and are one of the most promising places worth 
preempting for PV installation. 

First of all, previous construction characters of shopping centers facilitate the rooftop PV 
installation. According to literature [3], the availability rate of roof space is between 60% and 65% 
in shopping malls, but just 22% and 50% in residential buildings. That is to say, double or triple PV 
capacity can be installed in the former roofs, compared to the latter one. Secondly, the shopping 
malls need a large amount of electricity consumption daily, the most generation can feed the 
themselves-consumption. Besides, it is usually cited in downtown areas and populated areas. On the 
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one hand, there have a so complete transmission and distribution network that the surplus generation 
can efficiently feed local electricity consumption. On the other hand, brand advertisement and 
promotion of the PV manufacturer and the rising level of awareness of the citizen towards sustainable 
efforts can be greatly obtained for the large visitors there. Last but not least, facing the non-manageable 
feature through electricity generation in PV facilities, technical management measures need to be taken 
for protecting distribution system [4]. Fortunately, the shopping malls equipped professional staff and 
equipment for energy supply, security, air conditioning, and so on. So, the previous staff in shopping 
malls can condemn the whole new challenge to reduce the extra expenditure in PV management. 

It is desirable to adopt a proper methodology for evaluating the shopping mall PV plan in order 
to demonstrate the optimal possible selections for an investment. Since the following three problems 
existed in this issue, which are the interaction problems between the criteria, the loss of evaluation 
information in the conversion process, and the compensation problems between diverse criteria. 
There is no doubt that traditional multi-criteria decision-making (MCDM) method failed to solve the 
three problems at the same time. Nevertheless, the decision-making trial and evaluation laboratory 
(DEMATEL) method can help analyze the internal influence and connection behind each criterion, 
and the interval-valued neutrosophic set is accomplished in expressing the imperfect knowledge of 
experts group. Besides, an extended ELECTRE III method as an outstanding outranking method can 
succeed in avoiding the compensation problem, and the integrated method used maintained symmetry 
in the solar PV investment. Therefore, in this context, the integrated DEMATEL method and extended 
ELECTRE III method under interval-valued neutrosophic set environment for searching the optimal 
shopping mall solar PV plan has been devised. 


2. Literature Review 


MCDM (multi-criteria decision-making) has been successfully applied in energy planning 
problems. For example, a review of MCDM methods towards renewable energy development identified 
MCDM methods as one of the most suitable tools to finding optimal results concerned with energy 
planning progress in complex scenarios, including various indicators, and conflicting objectives and 
criteria [5]. For example, Fausto Cavallaro et al. use an intuitionistic fuzzy multi-criteria approach 
combined with fuzzy entropy to rank different solar-hybrid power plants successfully [6]. 

In a real case, it is different for decision makers to express preferences when facing inaccurate, 
uncertain, or incomplete information. Although the fuzz set, intuitionistic fuzzy sets, interval-valued 
intuitionistic fuzzy sets, and hesitant fuzzy sets can address the situation. However, when being 
asked the evaluation on a certain statement, the experts can use the interval-valued neutrosophic set 
(IVNNS) expressing the probability that the statement is true, false, and the degree of uncertainty can 
be accurately described, respectively [7]. The IVNNS, combined with outrank methods, has addressed 
many MCDM problems successfully [8]. For example, the IVNNS combined with VIKOR was applied 
to solve selection of location for a logistic terminal problem [9]. Hong-yu Zhang et al. developed 
two interval neutrosophic number aggregation operators and applied them to explore multi-criteria 
decision-making problems [10]. 

The independence of criteria remains in most of the MCDM methodologies. In recent years, lots of 
methods appeared to solve the problem, and the DEMATEL method is popularly used. According to 
the statistic censused in article [11], of the use of MCDM methods in hybrid MCDM methods, the top 
five methods are Analytical Network Process (ANP), DEMANTEL, Analytic Hierarchy Process(AHP), 
TOPSIS, and VIKOR. 

The DEMATEL methodology has been acknowledged as a proper tool for drawing the 
relationships concerning interdependencies and the intensity of interdependence between complex 
criteria in an evaluation index system [12,13]. As a powerful tool to describe the effect relationship, 
it help evaluate the enablers in solar power developments [11] and evaluate factors which influencing 
industries’ electric consumption [14]. The application of DEMATEL contributed to determining the 
weight coefficient of the evaluation criteria, and successfully helped identify the suitable locations 


Symmetry 2018, 10, 150 4 of 20 


for installation of wind farms. It is the DEMATEL method that helped the investors improve their 
decisions when there are many interrelated criteria. Thus, the connection relationship of the climate 
and economy criteria that exists in this study is of great need of the application of DEMATEL method. 

The commonly used outrank models are TOPSIS, AHP, ANP, preference ranking organization 
method for enrichment evaluations(PROMETHE), and ELECTRE, of which ELECTRE methods 
are preferred by decision makers in energy planning progress. Among the ELECTRE methods, 
ELECTRE III method conveys much more information than the ELECTRE I and ELECTRE II 
methods [15]. In the literature [16], economics of investment in the field of PV, an inclusive 
decision-making structure using ELECTRE III that would help photo voltaic (PV) system owners, 
bureaucrats, and the business communities to decide on PV technologies, financial support systems 
and business strategies were featured. ELECTRE III was used to structure a multi-criteria framework 
to evaluate the impact of different financial support policies on their attractiveness for domestic PV 
system deployment on a multinational level [17]. Due to the compensation problem in information 
processing, incomplete utilization of decision information, and information loss, ELECTRE III was 
chosen to build a framework for offshore wind farm site selection decision in the intuitionistic fuzzy 
environment. These literature studies improve the application of ELECTRE III method in the energy 
planning process, and terrify the effectiveness of evaluation in decision making progress [18]. 

In conclusion, based on the mentioned evolvement, the shopping mall PV plan evaluation result 
will be more scientific and reasonable than before. 


3. Decision Framework of SMPV Plan Selection 


The evaluation criteria are basic to the entire evaluation, so that they are of great importance to the 
shopping mall photovoltaic plan selection. In view of the special characteristics of photovoltaic 
plan and the shopping malls, six factors were taken into consideration, namely architectural 
elements, climate, photovoltaic array, economy, risk, contribution. Table 1 shows six criteria and 
twenty-one subcriteria. 


Table 1. Analysis of evaluation attributes of shopping centers photovoltaic plan selection. 


Criteria Subcriteria Resources 
al Roof pitch and orientation [2] 
(a) Architectural elements a2 Covering ratio [2] 
a3 PV roof space [2] 
bl Total investment [19] 
b2 Total profit [19] 
Pena: b3 Annual rate of return [20] 
b4 Payback year [20] 
cl Annual average solar radiation (kwh/ m?/ year) [19,21] 
(c) Climate c2 Land surface temperature (°C) [19,21] 
c3 Annual sunshine utilization hours (h) [19,21] 
dl Suitability of the local solar regime [22] 
: d2 PV area [2] 
d) Phot 
Ky Ep Hy aleeney d3 PV generation (yearly electricity generation) MWh/year [2,19] 
d4 Repair and clean rate [20] 
el Increase in local economy and employment [22] 
(e) Contribution e2 Publicity effects Own 
e3 Environment protection [23] 
fl Grid connection risk [19] 
(f) Risk 2 Rooftop ownership and occupancy disputes Own 
3 Bad climate [22] 


f4 Government subsidies reduction Own 
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3.1. Architectural Elements 


Not all the shopping malls are suitable to allocate the photovoltaics, and architectural elements 
are the primary intrinsic limitations. Steep roof pitch and wrong orientation will increase the difficulty 
of allocation and maintenance. In addition, building obstructions, and vegetation shading part of 
the space are not able to allow for the allocation of PV equipment, which can be measured by the 
covering ratio estimated by Equation (2). Last but not least, the photovoltaic roof space needed to 
be calculated by the Equation (1). Colmenar-Santos, Antonio et al. [15] assessed the photovoltaic 
potential in shopping malls by calculating the photovoltaic roof space. We decided to refer to this 
research method. 

RS PV = RS xa (1) 


RSpy is PV roof space, RS stand for roof space, «# is availability ratio 


Ragan = tan ¢ tan (as) 


: 2 
sin g tan (a;)" +sing tan gsin(75)" @) 

CR is covering ratio. («s)* solar altitude angle (7;)* g solar azimuth angle. The value is at nine 
o’clock on the winter solstice in each location. 


3.2. Climate 


Not all the locations of the shopping malls have the optimal climate for solar power generation. 
It is undoubted that the solar resource depends on local climate. Thus, the annual average solar 
radiation, land surface temperature, and annual sunshine utilization hours, are the four typical criteria 
to judge whether the local solar energy resource is in abundance. 


3.3. Photovoltaic Array 


It is well known that the performance of photovoltaic arrays will impact the electricity generation 
reliability and stability. The dust in the photovoltaic cell panel will reduce the solar energy conversion 
efficiency, meanwhile, since the photovoltaic cell panel damage and faults are directly related to the 
electricity supply reliability, the repair and clean rate should be pondered. In addition, it is worth 
concerning whether the specified type of the photovoltaic panel is absolutely suitable to the local 
solar regime. The total PV area affects the electricity generated which is calculated by Equation (3). 
PV generation (yearly electricity generation) is estimated by Equation (4). 


Apanel = RSpy x CR (3) 


E =H x Apanel X € X K (4) 


His the total yearly solar irradiation. A pane is the total area of PV panel. ¢ is the efficiency of the 
panel. « is the comprehensive facility performance efficiency, which is 0.8 [19], and 0.28 is the empirical 
conversion coefficient of the PV module area to the horizontal area. 


3.4. Economy 


It is beyond doubt that the economy of the shopping mall photovoltaic plans ought to be taken 
into account by the decision makers. There are plenty of studies to assess the financial aspects of 
the photovoltaic projects. Indrajit Das et al. [23] presented an investor-oriented planning model for 
optimum selection of solar PV investment decisions. Rodrigues Sandy et al. [24] conducted economic 
analysis of photovoltaic systems under China’s new regulation. The economic assessment methods 
there are so suitable and scientific that they are worth referring to in this paper. The significant 
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economic attributes we considered are pay pack period, total investment, total profit, and annual rate 
of return, which are calculated by Equations (5)-(7). 


Ppanel 


(5) 


I = C Apanel X 
panel 


I is the total investment, C is the average cost of building per W roof PV projects. Ppane is the 
max power pin (W) under STC situation of the solar panel, apane is the area of per photovoltaic panels. 


B=Prexttc+Ps xX ts -I x trc —O X ttc (6) 


B is the total profit, Pe is the electricity price buying from the power supply company, Ps is the 
electricity price subsidy, t;c is the time of PV projects life cycle, ts is the time subsidy lasting, O is the 


cost for operation and maintenance. 
B 


ROI = ——— 
2 Ix tre 


(7) 


ROT is the annual rate of return. ; 


> Pee SO 


(8) 


Tpg 
Tpg is the pay pack year. 


3.5. Contribution 


Although the environmental and social contributions the SMPV projects made may not be 
calculated explicitly as economic profit, there is no denying that these benefits result in increase in local 
economy and employment, and environment protection and publicity effects are worth the focus of 
attention. Particularly, the shopping mall holds a great number of visitors. On the one hand, it obtains, 
easily, the brand advertisement and promotion when PV equipment of a particular company occupies 
the rooftop of a large commercial building. On the other hand, it is effective to raise the level of 
awareness of the citizen towards renewable energy and sustainable efforts. 


3.6. Risk 


Expect that for the above attributes, the risk faced cannot be neglected. First of all, the government 
subsides policy is likely to change, and the impartiality, sufficiency, stability, and constancy of the 
subsidy is unable to be ensured. Secondly, the generating capacity is influenced by the climate heavily, 
so the profits will reduce when facing consecutive rainy days. Thirdly, the rooftop usage needs the 
allowance from all the owners, however, the rooftop ownership and occupancy disputes are a very 
common risk. Last but not least, the connected photovoltaic grid is unable to bring any benefits to the 
grid enterprise because of the intermittent power output. Thus, how long the support to photovoltaic 
grid connected from the grid enterprise can exist is uncertain. 

All in all, the SCPV plan alternatives ought to be appraised from architectural elements, climate, 
photovoltaic array performance, economy, risk, contribution attributes. The unique custom-made 
framework of criteria and subcriteria is set up in view of the actual SMPV plans and national conditions. 


4. Research Methodology 


A decision framework of SMPV selection has been proposed in this section, and there are 
four phases in this framework, as shown in Figure 2. The research framework is described in the 
following subsections. 
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Literature Experts 
Form experts group suvery opinion 
Solar RESOUMSES 2 rooftops Identification of criteria 
and other information ones 
and sub-criteria 
Possible selections of DEMATEL 
method 
SMPV 
Alternative shopping Criteria and sub- Influential network 
malls criteria weights among criteria 


Actual data of nvestigation and 
index quantitative calculation Eq(1- 
criteria 8) 
Define the discordance 
index 


Define the degree of 
credibility of the 
outranking relation 


IVNS score of 
; j creat IVNS score of 
Ranking orders obtained lhe ss ae qualitive criteria 


Compare analysis with the 
ranking odder from other 
methods 


Identification o 
the parameter on Experts opinion 
Eq.(25) 


Figure 2. The flowchart of the research methodology. DEMATEL: decision-making trial and 
evaluation laboratory. 


4.1. Preliminary Knowledge in the Neutrosophic Set Environment 


Due to the existence of many uncertainties in real decision-making problems, such as 
indeterminate and inconsistent, the neutrosophic set (NS) is used in the MCDM method, and definition 
of NS is introduced in this section. 


Definition 1 [25]. Let X be a space of objects with a generic element in X denoted by x. A NS A in X is 
defined using three functions: truth-membership function T 4 (x), indeterminacy-membership function I, (x) 
and falsity-membership function F4(x). These functions are real standard or nonstandard subsets of |O~,1*[, 
that is, Ta(x) : X >]0~,1*[, I4(x) : X >]07,1*[ and Fy(x) : X 4]07,1*[. And that the sum of T4(x), 
I4(x) and F4(x) satisfies the condition 0~ < supT,4(x) + supl,(x) + supF4(x) < 3°. 


Since the non-standard unit interval ]0~,1*| is hard to apply in practice, and the degree of truth, 
falsity, and indeterminacy about a certain statement could not be described precisely in the practical 
evaluation, the interval-valued neutrosophic set (IVNNS) of standard intervals has been proposed 
by Wang [26], and a few definitions and operations of IVNNS are introduced in the GPP technology 
selection MCDM problem. 
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Definition 2 [26]. Let X be a space of objects with a generic element in X denoted by x. An IVNNS A can be 
defined as 


A= {(x,Ta(x), Ia (x), Fa(x)) sx © X} (9) 


where Ta(x): X — [0,1], I4(x) : X > [0,1], Fa(x) : X — [0,1]. For each element x in X, these functions 
can be expressed as Ta(x) = [infT4(x),supT(x)] C€ [0,1], La(x) = [infla(x),suplI,(x)] C [0,1], Fa(x) = 
[infF.4(x),supF4(x)] C [0,1] and 0 < supT,(x) + supl,(x) + supF,4(x) < 3,x € X. For convenience, 
the interval-valued neutrosophic number (IVNN) can be expressed as @ = ((TE, T“], (12, I], [FE, FY]), 
and L, U represent the inferiors and superiors of IVNN respectively. 


ees oo L pu L yu L pu 5 — L 7u L yu L pu 
Definition 3 [10]. Leta = ([T4, 7H], [14,1], [FE EM]) and b = ([re, TU], [14,1], [PL FU] ) be two 
IVNNs, and A is a real number for not less than 0. Which operational rules can be expressed as follows 


woo =((e+ Th - ThTh TH + TY — Tere), ee], [etre eePU|) (10) 
i@b= bees ee [i +i- Ew + Pa), [Et + FL — pl.pl, PU 4 pu — Ful. pu) ) (11) 
aw = ((1-(1- TE)", - a — TH", [Cy ED”), [2] (12) 


w= (Te). (HF). [1- a - gy". -a- 8’), [1-G-Fey"1-0-)"]) a3) 


The original data of selection of SPPV are collected and processed in this phase. The alternative 
plans of the GPP project are evaluated by the experts, firstly according to the local technical condition 
data and practical experience. Then, the decision matrices are expressed in the form of IVNNs, 
which can handle incomplete and indeterminate information. Finally, a comprehensive decision 
matrix is formed based on interval-valued neutrosophic number weighted geometric operator 
(IVNNWG) operator. Let A; denote the technology alternatives (i = 1,2,---m), and Cj denote the 
criteria (j = 1,2,---n). It is assumed that a can be used to represent the evaluation value of attribute 
of alternative from every expert Ex(k = 1,2,---h). 


Definition 4 [27]. Let Ak — (a) be the IVNN-decision matrix of the k-th DM, k = 1,2,---h, 
mxn 


Wy 
~k L ru L yu 
and i= (74 73, eee, 


: Fi Ex 
J y y 
such that 


ij ij 


| ; An IVNNWG operator is a mapping: IVNN" + IVNN, 


IVNNWG, (a, a2, ne a) = (7) . 


k=h Wk k=h Wk k=h Wr k=h Wk k=h Wk k=h Wk 
= (C4) Cy FPR 0-4)" 0-8)" O-)-O-W)") 
k=1 \ 4 ka. \ Si k=1 aj k=1 ai k=1 aij k=1 aij 


where W = (W4,W2,-°-, wy)? represents the weight vector of DMs, satisfying Lt tg = 1, wx € [0,1]. 


4.2. Phase I Identification of Alternative SMPV Plans 


At this stage, a group of experts consisting of several doctorate engineers will be constituted by 
the investor. All the experts possess abundant working experience in solar energy investment field, 
and are specialized in solar photovoltaic and power grid technologies. 

More than twenty famous influential large-scale shopping malls located in those cities with both 
abundant sunshine and general policy subsidies need to be collected, and based on that information, 
less than ten alternative plans roughly screened out, based on the plentitude of documents and 
investigation. After that, the investigation will be carried out by the experts group, and involve 
meeting the Development and Reform Commission, the Meteorological Bureau, and the local power 
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supply companies, in order to gather information about solar resources, city planning and construction, 
local solar subsidy policies, distributed solar power planning, economic assessment, and approved 
shopping mall rooftops for construction. Lastly, there will be less than five of the most potential 
alternative places presented for the next evaluation. 


4.3. Phase II Determination of the Weights of Criteria Based on DEMATEL Method 


The importance of the criteria on GPP technology selection is different, so the DEMATEL method 
is used to decide the weight of criteria in this phase. The direct and indirect causal relations among 
criteria are considered in the DEMATEL method, and the subjective judgment of DMs is also considered. 
The steps of determining the weights based on the DEMATEL method are shown as follows [28]: 


Step 1. Determine the influence factors in the system 


The influence factors of GPP technology selection system are determined based on expert opinions 
and literature reviews, which is called the criteria, as shown in Table 1. 


Step 2. Construct the direct-relation matrix among the criteria 


The direct-relation matrix X = (Xpq),,,.,, is constructed in stages, where xpq is used to represent 
the degree of direct influence of pth criterion on gth criterion which is evaluated by the experts, and n 
is the number of criteria. 


Step 3. Normalize the direct-relation matrix 


The direct-relation matrix X = (Xpq) nxn is normalized into Y = (Yp4) nxn" A normalization 
factor [29] s is applied in the normalized calculation, and the normalized direct-relation matrix Y is 
calculated by using Equations (1) and (2). 


Y=s-X (15) 


s = Min 1 1 16 
(= q=1 Xpq) , Maxi <q<n (Dia 7) ( ) 


Step 4. Calculate the comprehensive-relation matrix. 


The comprehensive-relation matrix T is obtained by using Equation (3). 
T= ava (17) 
A=1 


where T = (tpq) . 4,7 P/4 = 1,2,-++n, and tpg is used to represent the degree of total influence of pth 
criterion on qth criterion. I represents for the identity matrix. 


Step 5. Determine the influence relation among criteria 


The influence degree and influenced degree of the criteria is determined after obtaining the 
comprehensive-relation matrix T. The sum of the row and column values of matrix T can be obtained 
by the Equations (4) and (5). The sum of row values of T, denoted by D, which represents the overall 
influence of a given criterion on other criteria. The sum of column values of T, denoted by R, which 
implies the overall influence of other criteria on a given criterion. 


D= (4p) nx = ae tpg (18) 


n 
oe (Ta) 1500 7 ae t pq (19) 
The causal diagram is obtained based on the D + R and D — R values. The D + R value indicates 
the importance of indicator in the SMPV plan selection system, the greater the D + R value, the more 
important the corresponding indicator is. On the other hand, the D — R value indicates the influence 
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between a certain indicator and the other indicators, which can be separated into cause and effect 
groups. The indicator, which has positive values of D — R, belongs to the cause group, and dispatches 
effects to the other indicators. Otherwise, the indicator, which has negative values of D — R, belongs to 
the effect group, and receives effects to the other indicators. 


Step 6. Determine the weight of criteria 


The criteria are represented by j, and satisfying j = 1,--- p,---q,---n, then the weights of criteria 
are determined based on the following Equations (6) and (7) [30]. 


wi! = [(dj +1))° + (dj —79)"| (20) 
w;! 
0) = Sr (1) 


where w rg denotes the relative importance of the indicators, and w j denotes the weights of the indicators 
in SMPV technology selection. 


4.4, Phase III Calculation IVNNs Performance Score 


ELECTRE is a family of methods used for choosing, sorting and ranking, multi-criteria problems. 
ELECTRE III was developed by Roy in 1978m, which is valued outranking relation. 

A = {a1,42,...,dn} is the finite set of alternatives, Y = {y1,¥2,...,Ym} is the finite set of criteria, 
y;(a) represents the performance of alternative a on criterion y; € Y. Assume that all the criteria 
are of the gain type, which means the greater the value, the better. q; is the indifference threshold, 
which represents two alternatives in terms of their evaluations on criterion y;. In general, qj is a 
function of attribute value q;(a;), which can be denoted as qj(y;(a;)); pj(y;(4i)) is preference threshold, 
which indicates that there is a clear strict preference of one alternative over the other in terms of their 
evaluations on criterion y;. In addition, v;(y;(a;)) is a veto threshold that indicates that the attribute 
value y;(a;) of scheme q; is lower than the attribute value y;(a,) of scheme a, and when it reaches or 
exceeds 0;(y;(a;)), it is not recognized that the a; is preferred to the ax. yj(a;) — yj(ax) which indicates 
the situation of preference of a; over a; for criterion C;. A weight w; expresses the relative importance 
of criterion y;, as it can be interpreted as the voting power of each criterion to the outranking relation. 

Where yj(a;) and yj(a,) are expressed in the form of IVNNs in the paper, that is, 
yas) = (TE, TH), (1, 2], (FE, FE]), y(n) = (TE, TE), (IE, I], [FE,EL]). In the calculation of 
Equation (8), let 


yj(ai) = (Th +7) — ip +1) — (FP + FH) (22) 
yj(ax) = (Te + TZ) — Ue + 2) — (Fe + FP) (23) 

and so 
yay) =e ey ee a ee eo) (24) 


bigmix (i € Op) 
Tj = Aj - (i E 0c, oo x 0) 
Ai(1 — ame) (i E Oc,amin 4 0) 


a 


Qij 7s 

ma: E 6 

oe B) (25) 
[ij = King (i E Ca x O)and 


Gij 


Hi(1 — gmax)(i € Oc, ae" # 0) 


Ej 
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The 6, A refer to the certainty parameters of the benefit criteria and cost criteria, respectively [23], 
while the ¢, » stand for the uncertainty parameters of the benefit criteria and cost criteria respectively, 
and they obey rule 0 < 6; + ¢; < 1,0<Aj+p; <1. 


4.5. Phase IV Calculation of Outranking Relation of IVNNs Based on Extended ELECTRE-III 


Step 1. Define the concordance index 


The concordance index c(aj,4,) that measures the strength of the coalition of criteria the support 
the hypothesis “is at least as good as”, c(a;,a,) is computed for each ordered pair a;,a, € A as follows: 


n 
C(aj, Ax) = hg (ai,4x)/ YW; (26) 
i=1 


and the partial concordance index c(a;,a;) is defined as 


0 (yj(4i) — yj(ax) < ajlyj(ai)]) 
ci(aj,a%) = 4 1 (yj(ai) — yj(ar) < ajlyj(ai)]) (27) 
yj (aj) —¥j (ax) —4jlyj (ai) (others) 
Ply (ai)]—9jly; (ai) 


Step 2. Define the discordance index 


The discordance index dj(a;,a,) is defined as follows: 


0 (yj(ax) — yj(ai) < —aylyj(ai 

dj(a;,ax)=4 1 (yj (ae) — yj(ai) = ojly;(ai)]) (28) 
yj (4x) —¥j (ai) +4; [yj (4) (others) 

oily; (ai) +4jly;(ai)] 


Step 3. Define the degree of credibility of the outranking relation 


The overall concordance and partial discordance indices are combined to obtain a valued 
outranking relation with credibility s(a;,a,) € [0,1] defined by: 


c(aji,a) TI Tae (others) (29) 
JES (ai4;) 


(aj, ar) (Vi, dj(aj, ax) < c(aj,a)) 
s(4aj, a) 


where j(a;,4,) is the set of criteria for which dj(a;,a,) > c(4j, ax). 

Step 4. Define the ranking of the alternatives 

bs(a; > a.) means the sum degree of credibility that alternative a; outranks all the other 
alternatives, and )'s(a, > a;) means the sum degree of credibility that all the other alternatives 


outrank alternative a;. Thus, AS(a;) represents the ranking of the alternative a;, and the higher the 
AS(a;) value is, the more superior the outranking order is. 


AS(a;) = > 8(a; > ag) — ) 8 (ay > aj) (k =1,2...n) (30) 
5. A Real Case Study 


A Chinese renewable energy investment company wants to build a shopping center rooftop 
photovoltaic power project. In order to seek the optimal shopping mall for rooftop photovoltaic power 
plants, furthermore, one must judge the weight and the influence network of the criteria. A group 
of experts consisting of three doctorate engineers (referred to as E1, E2, E3) was constituted by the 
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company. All three experts possess more than 15 years’ working experience in solar energy investment 
field, and are specialized in solar photovoltaic and power grid technologies. The collaboration of 
the all the experts was needed, thus, a pseudo-delphi method was applied in which each expert has 
no interaction. 

Considering the development target and the investment capacity of the company, the famous 
influential large-scale shopping malls located in those cities, with both abundant sunshine and 
general policy subsidies, have been roughly screened out based on the plentitude of documents 
and investigation. The investigation involved several potential shopping malls located in Beijing, 
Shanghai, Guangzhou, Chengdu, Hangzhou, and Nanjing. The experts group met the Development 
and Reform Commission, the Meteorological Bureau, and the local power supply companies, in order 
to gather information about solar resources, city planning and construction, local solar subsidy policies, 
distributed solar power planning, economic assessment, and approved shopping mall rooftops for 
construction. There are four potential shopping malls picked out after the first filter, and they are the 
Golden Resources shopping mall in Beijing, Super Brand Mall in Shanghai, Deji Plaza in Nanjing, 
Jiangsu province, The Mixc shopping mall in Shenzhen, Guangdong province (hereafter referred to as 
X1, X2, X3, X4), as shown in Figure 3. 


The alternative shopping malls geography distribution 


Golden Resources Shopping Mall 


Super Brand Mall 
Deji Plaza 


eT | 


The Mixe Shopping Mall a 


Figure 3. The alternative shopping malls geography distribution. 


Firstly, based on the evaluation criteria, the influence of each criteria to the other one criteria 
was accessed by the experts. Then, the three experts discussed with each other and obtained a 
consensus about the influence of each criteria, as shown in Table 2. According to the DEMANTEL 
method, the weight of criteria and subcriteria was calculated based on Equations (15)-(21), and shown 
in Table 2. For the intuitive and simple understanding and analysis of the criteria and subcriteria, 
Figures 4 and 5 were drawn. As shown in Figure 4, the horizontal axis represents the importance of 
a criteria, while the vertical axis indicates the influence between the criteria, and the arrow is from 
the sender of this influence to the receiver. As we can see, the (b) economy (0.286) obtained the most 
importance, but was vulnerable to other criteria. The (a) architectural element and (c) climate (0.88) 
seemed not particularly important, however, they had significant direct impacts to the other four 
criteria. The (f) risk (0.188), (d) photovoltaic array (0.162), and (e) contribution (0.151) were considered 
of medium importance, and among them, (f) and (d) had more of an impact, while (e) received more 
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impact. Horizontal histogram clearly and intuitively shows the weight of each subcriteria in Figure 5. 
It is obvious that the (b1) Total investment, (b2) Total profit, and (b3). Annual rate of return acquired the 
highest weight, in addition to the (f4) Government subsidies reduction, (e2) Publicity effects and (d2). 
PV area was considered to be less but also very important. Therefore, it can be imagined that the SMPV 
plan alternatives which obtained high scores in these criteria are more likely to win the competition. 


Gove 


LCI UNS 


ee 
Rooftop Ownership and o 


Environment protectin 
Increase in local economy and employment 
PV generation 

Suitability of the local solar regime 

Land surface temperature 


Pay back year — 
Total profit 


im ' J 


aaabbbbcccddddeeeffff 
123123412312341231234 


Roof pitch and G 


Figure 4. The weights of subcriteria. 


> D-R (c) Climate (0.125) 


elements 
(0.088) 


Figure 5. Influential network relationship map within systems. 


Secondly, there are two types of criteria, one is quantitative, and the other is qualitative. On one 
hand, for the quantitative subcriteria, searching from the NASA atmospheric science data center, 
the data for cl, c2, c3 were obtained. Through using the Google earth map, the roof space data 
were obtained. Then, according to Equations (1)—-(8), the data for a3, b1-b4, d2, d3 were estimated. 
Because the Equations (1)-(8) are just for rough estimate, these data were not highly accurate. 
Considering the uncertainty and fuzziness of the data, Equation (25) was used to turn the numerical 
value into an IVNN value. The performance scores of the quantitative subcriteria are shown in Table 3. 
On the other hand, for the qualitative subcriteria, the experts group devoted their efforts to investigate 
the alternative plans and evaluate the performance score for the subcriteria al, a2, d1, d4, el—e3, f1—f4. 
The performance scores of the qualitative subcriteria were shown in Table 4. In addition, the subcriteria 
were divided into positive and negative. The score of positive criteria higher and negative criteria 
lower means the alternative better. In this paper, subcriteria b1, d4, f1, f3, f4 are negative and the others 
are positive. 


14 of 20 


Symmetry 2018, 10, 150 


Table 2. The score of influence among each subcriteria. 
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b2 
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b3 


b4 
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c3 
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Table 3. IVNN performance scores of alternative SMPV plans on the quantitative subcriteria. 
Subcriteria S1 $2 S3 S4 
(a3) ([0.90,0.90],[0.10,0.10],[0.00,0.00]) — ({0.32,0.32],[0.04,0.04],[0.65,0.65]) — ([0.26,0.26],[0.03,0.03],[0.71,0.71])  ([0.42,0.42],[0.05,0.05],[0.54,0.54]) 
(b1) ([0.41,0.41],[0.10,0.10],[0.49,0.49])  ({0.70,0.70],[0.08,0.06],[0.23,0.25])  ([0.26,0.26],[0.03,0.03],[0.71,0.71])  ([0.50,0.50],[0.06,0.06],[0.44,0.44]) 
(b2) ([0.90,0.90],[0.03,0.03],[0.07,0.07])  ({0.40,0.40],[0.04,0.06],[0.56,0.54]) — ([0.90,0.90],[0.10,0.10],[0.00,0.00]) — ([0.53,0.53],[0.05,0.05],[0.42,0.42]) 
(b3) ([0.90,0.90],[0.03,0.03],[0.07,0.07])  ({0.65,0.65],[0.07,0.07],[0.28,0.28])  ([0.53,0.53],[0.06,0.06],[0.41,0.41])  ([0.63,0.63],[0.07,0.07],[0.30,0.30]) 
(b4) ([0.90,0.90],[0.03,0.03],[0.07,0.07])  ({0.72,0.72],[0.08,0.08],[0.20,0.20]) — ([0.60,0.60],[0.07,0.07],[0.33,0.33])  ([0.60,0.60],[0.07,0.07],[0.33,0.33]) 
(cl) ([0.90,0.90],[0.03,0.03],[0.07,0.07])  ({0.79,0.79],[0.10,0.10],[0.11,0.11])  ({0.81,0.81],[0.09,0.09],[0.10,0.10]) — ([0.81,0.81],[0.09,0.09],[0.09,0.09]) 
(c2) ([0.53,0.53],[0.06,0.06],[0.41,0.41])  ({0.70,0.70],[0.08,0.08],[0.22,0.22])  ([0.62,0.62],[0.07,0.07],[0.31,0.31])  ([0.90,0.90],[0.10,0.10],[0.00,0.00]) 
(c3) ([0.90,0.90],[0.03,0.03],[0.07,0.07])  ({0.71,0.71],[0.08,0.08],[0.21,0.21])  ([0.69,0.69],[0.08,0.08],[0.23,0.23])  ([0.73,0.73],[0.08,0.08],[0.19,0.19]) 
(d2) ([0.90,0.90],[0.03,0.03],[0.07,0.07])  ({0.53,0.53],[0.06,0.06],[0.41,0.41])  ([0.41,0.41],[0.05,0.05],[0.54,0.54]) — ([0.74,0.74],[0.08,0.08],[0.17,0.17]) 
(d3) ([0.90,0.90],[0.03,0.03],[0.07,0.07])  ({0.47,0.47],[0.05,0.05],[0.48,0.48])  ([0.37,0.37],[0.04,0.04],[0.59,0.59]) — ([0.67,0.67],[0.07,0.07],[0.25,0.25]) 
Table 4. IVNN performance scores of alternative SMPV plans on the qualitative subcriteria. 
Subcriteria X1 X2 X3 X4 
(al) ([0.45,0.56],[0.18,0.30],[0.13,0.50]) — ({0.68,0.79],[0.18,0.30],[0.13,0.50]) — ([0.56,0.68],[0.18,0.24],[0.38,0.67]) — ([0.79,0.90],[0.12,0.18],[0.25,0.50]) 
(a2) ([0.50,0.50],[0.06,0.06],[0.44,0.44])  ({0.85,0.85],[0.09,0.09],[0.06,0.06])  ({0.80,0.80],[0.09,0.09],[0.12,0.12])  ([0.90,0.90],[0.10,0.10],[0.00,0.00]) 
(d1) ([0.80,0.90],[0.15,0.23],[0.13,0.33])  ({0.50,0.60],[0.23,0.30],[0.38,0.67])  ([0.40,0.60],[0.15,0.23],[0.25,0.50]) — ([0.60,0.80],[0.23,0.30],[0.25,0.50]) 
(d4) ([0.20,0.23],[0.12,0.20],[0.25,0.33])  ({0.45,0.60],[0.20,0.30],[0.13,0.17])  ({0.60,0.90],[0.15,0.20],[0.13,0.14])  ([0.30,0.45],[0.10,0.12],[0.20,0.50]) 
(e1) ([0.70,0.90],[0.08,0.15],[0.21,0.33])  ({0.50,0.60],[0.23,0.30],[0.36,0.58]) — ({0.30,0.50],[0.08,0.15],[0.43,0.58]) — ([0.60,0.80],[0.15,0.30],[0.14,0.25]) 
(e2) ([0.50,0.60],[0.23,0.30],[0.33,0.75]) — ({0.70,0.80],[0.15,0.23],[0.33,0.75]) — ({0.70,0.90],[0.15,0.23],[0.17,0.50]) — ([0.80,0.90],[0.15,0.23],[0.17,0.50]) 
(e3) ([0.79,0.90],[0.12,0.18],[0.17,0.50]) — ({0.68,0.90],[0.06,0.12],[0.17,0.75]) — ([0.68,0.79],[0.18,0.24],[0.33,0.75]) — ([0.68,0.79],[0.24,0.30],[0.33,0.75]) 
(f1) ([0.30,0.45],[0.20,0.26],[0.33,0.43])  ({0.45,0.90],[0.26,0.30],[0.38,0.43])  ([0.23,0.45],[0.20,0.23],[0.38,0.43])  ([0.30,0.90],[0.23,0.26],[0.38,0.50]) 
(f2) ([0.45,0.68],[0.08,0.15],[0.44,0.56])  ({0.45,0.68],[0.15,0.23],[0.39,0.56]) — ([0.23,0.45],[0.15,0.23],[0.44,0.56])  ([0.68,0.90],[0.23,0.30],[0.33,0.44]) 
(£3) ([0.23,0.26],[0.12,0.20],[0.25,0.33])  ({0.60,0.90],[0.20,0.30],[0.13,0.17])  ([0.60,0.90],[0.15,0.20],[0.13,0.14])  ([0.30,0.45],[0.10,0.12],[0.20,0.50]) 
(f4) ([0.30,0.34],[0.15,0.30],[0.25,0.50]) — ({0.30,0.39],[0.15,0.30],[0.17,0.25])  ([0.68,0.90],[0.03,0.04],[0.13,0.25])  ([0.68,0.90],[0.04,0.05],[0.17,0.25]) 
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Thirdly, based on the weight of the sub-criteria and the IVNN scores of each alternative on each 
subcriteria, the final composite scores were calculated by improved ELECTRIC III method. After being 
told that q; is the indifference threshold, p; is the preference threshold and v; is the veto threshold, 
the experts group suggested that )°s(a;,a,) — Vs(ag,a;) (k =1,2...n), respectively. The concordance 
index c(X;, X;) and the partial concordance index c;(X;, X;) were calculated by Equations (26) and (27), 
as shown in Table 5. The discordance index dj(X;, Xx) was achieved by Equation (28), as shown in 
Table 6. Then, overall concordance and partial discordance indices was obtained by Equation (29) 
as shown in Table 7. Finally, the degree of credibility of the outranking relation was calculated by 
Equation (30), and the rankings of alternative X1, X2, X3, X4 was shown in Table 8. 

From Table 8, the SMPV plan X1 of the Golden Resources shopping mall in Beijing is the optimal 
selection. The alternative X1 is particularly superior to other alternative plans in terms of the economy, 
photovoltaic array, and contribution criteria, while these three criteria weighed more than a half of the 
entire criteria weights, so there is no doubt that plan X1 obtained the best position. However, plan X1 
performs badly in the risk and architectural elements criteria. Respectively, Plan X2 have strength on 
the economy, but are weak on photovoltaic criteria. Yet, plan X2 is much better than plan X3 and X4, 
so that it can be the stand-by choice. 


Table 5. The concordance index and the partial concordance index for each pair of SMPV plans. 


al a2 a3 b1 b2 b3 b4 cl c2 c3 di 


c(X1>X2) 0.00 0.00 1.00 0.00 1.00 062 043 0.24 0.00 046 0.83 
c(X1>X3) 0.05 0.00 1.00 0.36 0.00 0.92 0.74 0.20 0.00 0.50 0.61 
c(X1>X4) 0.00 0.00 1.00 0.00 0.93 0.68 0.74 0.19 0.00 041 0.45 
c(X2>X1) 0.26 0.85 0.00 0.70 0.00 0.00 0.00 0.00 041 0.00 0.00 
c(X2>X2) 034 010 O11 1.00 0.00 0.26 0.28 0.00 0.19 0.01 0.00 
c(X2 > X3) 0.00 0.00 0.00 047 0.00 0.02 0.28 0.00 0.00 0.00 0.00 
c(X3>X1) 0.00 0.72 0.00 0.00 0.00 0.00 0.00 0.00 0.20 0.00 0.00 
c(X3 > X2) 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.01 0.00 0.00 0.19 
c(X3 > X4) 0.00 0.00 0.00 0.00 0.93 0.00 0.00 0.00 0.00 0.00 0.00 
c(X4>X1) 044 0.99 0.00 0.19 0.00 0.00 0.00 0.00 0.92 0.00 0.00 
c(X4>X2) 015 O11 0.22 0.00 0.30 0.00 0.00 0.02 048 0.02 0.35 
c(X4>X3) 0.52 0.24 037 059 0.00 0.21 0.00 0.00 069 0.05 0.13 


d2 d3 d4 el e2 e3 f1 f2 £3 £4 C 


c(X1>X2) 0.91 1.00 0.00 0.74 0.00 0.12 0.00 0.03 0.00 0.00 0.38 
c(X1>X3) 1.00 1.00 0.00 0.78 0.00 046 0.02 035 0.00 0.00 0.40 
c(X1 > X4) 037 055 0.00 014 0.00 054 0.00 0.00 0.00 0.00 0.32 
c(X2>X1) 0.00 0.00 044 0.00 0.32 0.00 0.26 0.00 0.70 0.21 0.18 
c(X2>X2) 0.28 0.23 0.00 0.01 0.00 030 032 0.29 0.00 0.00 0.19 
c(X2>X3) 0.00 0.00 0.24 0.00 0.00 038 0.06 0.00 0.53 0.00 0.11 
c(X3>X1) 0.00 0.00 085 0.00 0.66 0.00 0.00 0.00 081 1.00 0.19 
c(X3 > X2) 0.00 0.00 0.37 0.00 0.30 0.00 0.00 0.00 0.08 0.81 0.16 
c(X3 > X4) 0.00 0.00 0.65 0.00 0.00 0.05 0.00 0.00 065 0.00 0.13 
c(X4>X1) 0.00 0.00 0.17 0.00 0.72 0.00 017 0.22 0.13 1.00 0.20 


c(X4>X2) 051 049 0.00 0.56 0.37 0.00 0.00 0.28 0.00 0.78 0.21 
c(X4> X3) 0.82 0.75 0.00 0.61 0.03 0.00 0.22 060 0.00 0.00 0.25 
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Table 6. The discordance index for each pair of SMPV plans. 

d(Xi > Xj) all a2 a3 b1 b2 b3 b4 c1 c2 c3 di 
d(X1>xX2) 0.00 0.00 1.00 000 091 047 034 0.21 0.00 0.36 0.61 
d(X1>xX3) 0.08 0.00 1.00 030 0.02 068 056 019 0.00 0.39 0.46 
d(X1>xX4) 0.00 0.00 088 0.00 068 051 056 O17 0.00 0.33 0.35 
d(X2>X1) 0.22 063 0.00 052 0.00 0.00 0.00 0.00 0.33 0.00 0.00 
d(X2>X2) 0.28 0.11 0.12 0.80 0.00 0.23 O24 0.00 0.17 20.05 0.00 
d(X2 > X3) 0.00 0.00 0.00 0.37 0.00 0.06 0.24 0.00 0.00 0.00 0.00 
d(X3>X1) 0.00 054 0.00 0.00 0.02 0.00 0.00 0.00 0.18 0.00 0.00 
d(X3>X2) 0.00 0.00 0.00 0.00 0.91 0.00 0.00 0.05 0.00 0.00 0.17 
d(xX3 > X4) 0.00 0.00 0.00 0.00 068 0.00 0.02 0.01 0.00 0.00 0.00 
d(X4>X1) 035 0.72 0.00 018 0.00 0.00 0.00 0.00 068 0.00 0.00 
d(xX4>X2) 015 012 0.20 0.00 025 0.00 0.00 0.06 037 0.06 0.29 
d(xX4>xX3) 041 021 030 045 0.00 019 002 0.03 052 0.08 0.13 
d(Xi > Xj) d2 d3 d4 el e2 e3 f1 f2 £3 f4 
d(X1>X2) 0.67 0.78 0.00 055 0.00 0.13 0.000 0.06 0.00 0.00 
d(X1>xX3) 089 0.96 0.00 058 0.00 0.36 0.06 0.29 0.00 0.00 
d(X1>X4) 030 043 0.00 014 0.00 041 0.00 0.00 0.00 0.00 
d(X2>X1) 0.00 0.00 035 0.00 0.27 0.00 0.23 0.00 052 0.19 
d(X2>X2) 0.24 0.20 0.00 0.05 0.00 0.25 0.26 0.25 0.00 0.00 

d(X2 > X3) 0.00 0.00 0.21 0.00 0.00 031 0.09 0.00 041 0.00 
d(x3>X1) 0.00 0.00 063 0.00 050 0.00 0.00 0.00 060 0.77 
d(X3>X2) 0.00 0.00 030 0.00 025 0.00 0.00 0.00 0.10 0.60 

d(xX3 > X4) 0.00 0.00 049 0.00 0.00 0.08 0.00 0.00 049 0.05 
d(x4>X1) 0.00 0.00 016 0.00 054 0.00 016 019 0.13 0.75 
d(X4>X2) 039 0.38 0.00 043 0.30 0.00 0.00 0.24 0.00 0.58 
d(x4>X3) 0.61 056 0.00 046 0.07 0.00 0.20 046 0.00 0.00 


Table 7. The overall concordance and partial discordance indices for each pair of SMPV plans. 


s(Xi > Xj) al a2 a3 b1 b2 b3 b4 c1 c2 c3 di 
s(X1>X2) 162 167 0.00 167 015 088 1.10 132 167 1.07 0.64 
s(X1>X3) 154 167 0.00 118 163 054 0.74 136 1.67 1.02 0.90 
s(X1>X4) 148 148 018 148 047 0.73 0.66 1.22 148 1.00 0.96 
s(XX2>X1) 0.95 045 1.22 058 1.22 1.22 1.22 1.22 082 1.22 1.22 
s(X2>X2) 088 1.09 1.08 025 1.23 0.95 094 1.23 102 1.17 1.23 
s(X2>X3) 112 112 112 O70 112 105 O85 112 112 4112 «1.12 
sX3>X1) 1.23 O57 1.23 123 1.20 1.23 1.23 123 101 1.23 = 1.23 
s(XX3>X2) 1.19 119 1.19 119 O10 119 119 113 119 119 0.98 
s(XX3>X4) 114 114 %114 +1114 O36 114 %112 #2113 #114 «-21214~«21214 
s(X4>X1) 082 034 1.25 103 1.25 1.25 1.25 125 040 1.25 1.25 
s(X4>X2) 1.09 1.12 1.02 127 095 1.27 1.27 1.20 080 1.20 0.91 
s(X4>X3) 0.79 1.06 094 073 133 1.08 130 1.29 064 1.22 1.16 
s(Xi > Xj) d2 d3 d4 el e2 e3 f1 f2 £3 £4 s 

s(X1>X2) 055 036 167 0.75 167 145 167 1.56 1.67 1.67 0.38 
s(X1>X3) 0.19 0.06 167 069 1.67 1.07 1.57 1.19 167 1.67 0.40 
s(X1>X4) 104 085 148 1.27 148 087 148 148 148 148 0.32 
s(XX2>X1) 1.22 1.22 0.79 1.22 089 1.22 0.94 1.22 0.58 0.98 0.18 
s(X2>X2) 0.94 0.98 1.23 116 1.23 0.92 091 0.93 1.23 1.23 0.19 
s(X2>X3) 1.12 112 088 112 112 O78 102 112 066 112 O11 
sX3>X1) 1.23 123 046 1.23 062 1.23 123 123 049 0.28 0.19 
s(XX3>X2) 1.19 1.19 083 1.19 089 1.19 119 119 1.07 047 0.16 
s(X3>X4) 114 114 O58 %114 +114 +106 114 +114 058 1.09 0.13 
sX4>X1) 1.25 1.25 105 125 058 1.25 1.05 1.01 1.08 0.32 0.20 
s(X4>X2) 0.77) 0.79 1.27 0.72 0.89 1.27 1.27 0.97 1.27 054 0.21 
s(X4>X3) 052 059 133 O71 1.24 133 1.07 072 133 1.33 0.25 
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Table 8. Final composite scores and rankings of alternative SMPV plans. 


Xi Ys(aj > ax) Yis(ax > aj) AS(aj) 
X1 0.73 0.29 0.44 
X2 0.47 0.37 0.1 
X38 032 0.39 —0.07 
X4 0.56 1.02 —0.47 


The ranking order using ELECTRI III is X1 > X2 > X3 > X4. In order to check the validity of the 
results, TOPSIS and VIKOR methods were used to reorder the alternative SMPV plans as shown in 
Table 9. The result obtained by TOPSIS method is X1 > X4 > X3 > X2, while the result achieved by 
VIKOR method is X1 > X4 > X3 > X2. from these three rankings, alternative plan X1 is the optimal 
selection, no matter what method was used. That is to say, the alternative X1 is much better than 
the remaining alternatives, and there is no doubt in choosing X1 first. However, the rankings for X2, 
X3, and X4 are different between these three methods. There is the veto threshold, which indicates 
when the value of alternative Xi is lower than the value of alternative Xj, and the lower value exceeds 
the veto threshold, and it is not recognized that Xi is preferred to the Xj in general. However, it is 
the other term in TOPSIS method and VIKOR method, and some really bad performance in a certain 
criterion can be tolerated and remedied by other good performances in other criteria. In that case, 
an alternative with some fatal defect in a certain criterion of an alternative may be neglected, which 
leads to an unsatisfactory selection. When an alternative is vetoed better than the other alternative in 
ELECTRE II, it can still come out in front in the TOPSIS and VIKOR methods. That why the X2, X3, 
X4 ranked differently in TOPSIS and VIKOR methods. 


Table 9. The rankings of SMPV plans using TOPSIS and VIKOR. 


Method TOPSIS VIKOR 
yt y- Cc Rankings s rt QO Rankings 
X1 0.62 1.43 0.70 1 —2.57 0.12 0.00 1 
X2 1.17 0.47 0.29 4 1.88 1.18 0.922 4 
X3 1.35 0.71 0.35 3 —174 1.06 0.526 2 
X4 0.76 0.74 0.49 2 2.74 1.01 0.919 3 


6. Conclusions 


The selection of SMPV plan is crucial to the entire life of SMPV project. Although there has been 
some research on this issue, several questions still need addressing. Firstly, the interaction of the criteria 
lay in the evaluation criteria. Secondly, the loss of evaluation information exited in the information 
conversion process. Thirdly, the compensation problem between best and worst performance in 
diverse criteria was not easily to avoided. 

In this paper, an integrated MCDM framework was proposed to address the SMPV plan selection 
problem. First of all, the compositive evaluation index was constructed, and the application of 
DEMATEL method helped analyze the internal influence and connection behind each criterion. 
From the influential network-relationship map, we discovered that the criteria (b) economy obtained 
the most importance but was vulnerable to other criteria as well as the (a) architectural element and (c) 
climate had significant direct impacts to the other four criteria. These three criteria should be the first for 
the decision maker to consider when selecting the SMPV plan. Then, the interval-valued neutrosophic 
set is utilized to express the imperfect knowledge of experts group. Since the application of IVNNS, 
the experts can clearly express their evaluation information, including their certainty, uncertainty, 
as well as hesitation attitude. Following this, an extended ELECTRE III method as an outstanding 
outranking method was applied, and it succeed in avoiding the compensation problem and obtaining 
the scientific result. In the case of China, the integrated method has been successfully applied to select 
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the SMPV plan X1 as the optimal selection which is particularly superior to other alternative plans in 
terms of the economy, photovoltaic array, and contribution criteria. Also, the integrated method used 
maintained symmetry in the solar PV investment. Last but not least, a comparative analysis using 
TOPSIS method and VIKOR method was carried out, and alternative plan X1 ranks first at the same. 
The outcome certified the correctness and rationality of the results obtained from this paper. 

Therefore, this study has not only served to evaluate the SMPV plans, it has also demonstrated 
how it is possible to combine IVNNS, DEMATEL method, and ELECTRE III method for application in 
handling MCDM problems in the field of solar energy. 
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